Optical phase estimation in the presence of phase-diffusion 
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The measurement problem for the optical phase has been traditionally attacked for noiseless schemes or in 
the presence of amplitude or detection noise. Here we address estimation of phase in the presence of phase 
diffusion and evaluate the ultimate quantum limits to precision for phase-shifted Gaussian states. We look for 
the optimal detection scheme and derive approximate scaling laws for the quantum Fisher information and the 
optimal squeezing fraction in terms of the total energy and the amount of noise. We also found that homodyne 
detection is a nearly optimal detection scheme in the limit of very small and large noise. 
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The estimation of the optical phase in quantum mechanics 
is a longstanding problem with both fundamental and tech- 
nological implications. The attempts to define a Hermitian 
phase-operator conjugate to the number operator are the sub- 
ject of an extensive literature ^ and several experimental pro- 
tocols to estimate the value of the optical phase have been 
also proposed and demonstrated, in particular using different 
quantum strategies and interferometric setups yHZl], which 
have been shown to beat the standard quantum limit and to 
attain the Heisenberg limit H - H | . More recently the ultimate 
bounds to precision of phase estimation with Gaussian states 
have been discussed in details lll2i 11311 using local quantum 
estimation theory. Squeezed vacuum state has been shown to 
be the most sensitive for a given average photon number and 
two adaptive local measurement schemes have been proposed 
to attain the Heisenberg limit asymptotically. 

The estimation of the optical phase, besides the fundamen- 
tal interest, is also relevant for optical communication scheme 
where information is encoded in the phase of laser pulses that 
must travel long distances between the sender and receiver. 
In such a context the receiver has to decode the information 
carried by the optical wave-packets which will be unavoid- 
ably degraded by different sources of noise, which have to 
be duly taken into account in the quantum estimation prob- 
lem. As a matter of fact, only amplitude and/or detection 
noise have been taken into account in the analysis of quantum 
phase estimation, e.g. imperfect photodetection in the mea- 
surement stage, or amplitude noise in interferometric setups 
13-0 ■ The role of phase-diffusive noise in phase-estimation 
have been investigated for qubit 1 19, i^l and in part for con- 



densate systems 112 IL 12211 . while no similar analysis have been 
done for a continuous variable system. Phase-diffusive noise 
is the most detrimental for phase-estimation since it destroys 
the off-diagonal elements of the density matrix. Moreover, 
any quantum state that is unaffected by phase-diffusion, is 
also invariant under a phase-shift, and thus is totally useless 
for phase-estimation. 

In this letter we address for the first time phase estimation in 
the presence of phase diffusion, seek for the optimal Gaussian 
states, and evaluate the ultimate quantum limits to precision 



of phase estimation. We also investigate whether the ultimate 
performances may be achieved with feasible detection scheme 
and found that homodyne detection is nearly optimal for very 
small and large amount of noise. 

Let us start by a pico review of local quantum estimation 



theory 023142611 . When a physical parameter is not directly ac- 
cessible one has to resort to indirect measurements. Let us 
denote by the quantity of interest, X the measured observ- 
able, and X = (^ii ■ ■ • : ^m) the observed sample. The es- 
timation problem amounts to find an estimator, that is a map 
4> — 4>{x) from the set of the outcomes to the space of param- 
eters. Classically, optimal estimators are those saturating the 
Cramer-Rao inequality Var(0) > [MF{4))]^^ which bounds 
from below the variance Var(0) ~ E[(fP'] — E[(f)Y of any un- 
biased estimator of the parameter 0. In the Cramer-Rao in- 
equality, M is the number of measurements and is the 
Fisher Information (FI) = j dxp{x\4i)\dij,\Tip{x\4i)\ 

where p(a;|0) is the conditional probability of obtaining the 
value X when the parameter has the value (f>. The quantum 
analogue of the Cramer-Rao bound is obtained starting from 
the Born rule p{x\(j)) = Trpj^p^] where {lix} is the proba- 
bility operator-valued measure (POVM) describing the mea- 
surement and the density operator, labeled by the param- 
eter of interest. Upon introducing the Symmetric Logarith- 
mic Derivative (SLD) as the operator satisfying 2d^g^ = 
L^g^ + QtjjL^ one proves that the FI is upper bounded by the 
Quantum Fisher Information (QFI) jH F{4>) < H{(j)) = 
Tr[p^L|]. In turn, the ultimate limit to precision is given 
by the quantum Cramer-Rao bound Var((/)) > [AIH{<f>)]^^. 
The family of states we are going to deal is a unitary one 
= U^goUl = EfeAfe I Afe ((/.)) (Afc ((/))!, where |Afe(0)) = 
Ucj,\Xk) and C/^ = exp{—i(j)G} describes a phase-shift with 
the single-mode number operator G = a^a as the generator 
In this case the SLD may be written as = J/^LqC/^, where 
Lq is independent on cj). The corresponding QFI does not de- 
pend on the parameter (j), and reads 



H = TvigoQ = 2 



(A,i — A,„)^ 
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Phase-diffusion for a continuous-variable system is de- 
scribed by the master equation g = r£[a^a]p where 
C[0]q = 20(?0t _ otop _ gOto. The solution for 
an initial state g(0) is given by g(t) = A/'A(g(0)) = 
En,,ne-^'("-'")'e„,™(0)|n)(m| where A = Ft, T is the 
noise amplitude and (?„,m(0) = (n|p(0)|TO). The diagonal el- 
ements Q are left unchanged and, in turn, energy is conserved, 
whereas the off-diagonal ones are progressively destroyed. 

We assume that phase noise occurs between the applica- 
tion of the phase-shift and the detection of the signal, and 
address quantum estimation of a phase-shift applied to pure 
single-mode Gaussian states j^/'c) = D{a)S{r)\Q) where 
S{r) ~ exp{(r/2)(a^ — a^^) is the squeezing operators, 
D{a) = exp{Q!(a^ ~ o)} the displacement operator, being 
r,a G R. The input state is firstly phase-shifted by apply- 
ing the unitary operator U^, where is the unknown phase- 
shift, and then, before being measured, it undergoes phase- 
diffusion. Our aim is to determine the ultimate bound to pre- 
cision for a generic pure Gaussian probe and then look for the 
optimal one by maximizing the QFI over the state parameters. 

The mixed non-Gaussian state that is being measured is 
given by 

Q4,{t)^NA{U^\^G){i'G\Ul) = U4,MA{\i^G){M)Ul, 

where the second equality holds since the superoperator 
£[a^a] and the phase-shift operator [/^ commute. Be- 
cause of this fact our estimation problem corresponds 
to the case of a unitary family described above, with 
the input mixed state given by A/'a(|V'g)(^'g|)- In or- 
der to evaluate the corresponding QFI one writes £»0 
in its diagonal form = Y.n^ri\Ki{4>)){K{4>)\ = 
Y.n'^riU^\\n){Xn\U'^^, whcrc |A„(0)) and |A„) are respec- 
tively the eigenvectors of g^f, and of A/'a(|V''g)('0g|) cor- 
responding to the eigenvalues A„, which are in fact left 
unchanged by the phase-shift operation. By decomposing 

An) = X^fe Tnklk) in the Fock basis and by substituting this 
into the eigenvalues equation MAiH'c) {''PgDI^ti) ~ A„|A„) 
wehecve {n\ijG){i^G\l)e^^^^"^''^^rgk = Xgrqu^n. More- 
over, since a^a|A„) = ^''"^ s*'^'^!^)' that: 

(A„i|a1^a|A„)p = Efc ^ ^mfc J'nfc After evaluating the QFI 
using the above formulas one sees that it depends only on the 
eigenvalues A„ and on the components of the eigenvectors r„fc 
which, being a unitary parameter, do not depend on the pa- 
rameter itself. The explicit values of A„ and r„fe have been 
obtained by performing numerical diagonalization. 

Upon inspecting the solution of the master equation one 
sees that the vanishing of the off-diagonal matrix elements is 
governed by the product between and the squared differ- 
ence between the Fock indices {n — m)^. Besides, for a pure 
Gaussian state, the presence of non-zero (non-negligible) off- 
diagonal elements is somehow ruled by the average photon- 
number N = (a^a) and thus we roughly expect the QFI to 
somehow depend on the quantity ^ = A^^A. Pure Gaussian 
states may conveniently parametrized by the average photon 
number of photons N and of the corresponding squeezing 



fraction /?, in formula A'' = sinh^ r+jap and /? = sinh^ r/N, 
and thus the QFI will be function of the three parameters N, f3 
and A. 

We start our analysis by evaluating the QFI at fixed noise 
A. We consider four values of the maximum energy A^max = 
(a^^a)niax = {10,15,20,30} (with 10 steps on intermediate 
energies N) and different values of the noise parameter A. 
The values of A are chosen such that we can find points corre- 
sponding to fixed values of The curves are built by looking 
for the optimal pure Gaussian state, i.e. maximizing the QFI 
as a function of the squeezing fraction /3, for any fixed value 
of the energy N and of the noise parameter A. 




N 



2 4 6 8 10 2 4 6 8 10 12 14 



^opt ^opt 




N 



5 10 15 20 5 10 15 20 25 30 



FIG. 1 : (Color online) Optimal squeezing fraction /3 as a function of 
the average photon number and for different values of A^. (Top 
left): from top to bottom = {4.5 x 10-^ 4.5 x 10"*, 4.5 x 10-^ 
4.5 X 10"^}. (Top right): from top to bottom = {2.0 x 10"^ 
2.0 x 10""*, 2.0 X 10"^, 2.0 X 10"^}. (Bottom left): from top 
to bottom = {1.125 x 10"^ 1.125 x 10"*, 1.125 x 10"^ 
1.125 x 10"^}. (Bottom right): from top to bottom A ^ = {5.0 x 
10"^ 5.0 x 10~^ 5.0 x 10"*, 5.0 x 10"^}. 

The values of the optimal squeezing fraction /?opt = 
/Jopt(Af, A) and of the coiTesponding QFI i7(iV, /^opt, A) 
have been numerically evaluated and are reported in Fig. [T] 
and Fig. |2] respectively. As we can see in Fig. [1] for a low 
level noise the squeezing fraction is almost equal to one. In 
particular, in each plot, for the lowest value of A, we ob- 
tain /3opt{N, A) = 1 and thus the optimal probe state is the 
squeezed vacuum state, as it happens in the noiseless case 
il2ll . As far as the noise A increases the squeezing frac- 
tion decreases as a function of the average number of pho- 
tons. This means that for increasing values of the noise and 
of the energy, it is more convenient to employ the energy in 
increasing the coherent amplitude rather than the squeezing 
of the probe. Let us now focus on the behavior of the QFI 
H{N, /?opt, A). In the left panel of Fig. |2l we report the typ- 
ical behavior of the QFI as a function of N and for different 
values of A. The QFI increases by increasing the average pho- 
ton number A^, and decreases with the noise parameter A. For 
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the lowest value of A, we also observe that the noiseless limit 
H{N,I3 = 1, A = 0) = 8{N^ + TV) m is approached, at 
least for N not too large. 



H (10^) 




A. We notice that the H{N, A) is decreasing exponentially 
with the phase noise and that higher values of N correspond 
to higher values of H. 
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FIG. 2: (Color online) Left panel: QFI of optimized pure input 
Gaus.sian states as a function of the average photon number N and 
for different values of the noise parameter A; from top to bottom 
= {5.0 X 10"^ 5.0 X 10"^ 5.0 x lO"'*}. The black dot- 
ted line is the QFI for the noiseless case H{N, P = 1, A = 0) — 
8{N'^ + N). Right panel: (orange points) — ln7(^) as a function of 
InC, C = NA, with iV < 30 and 10"^ < A < 1. The black curve 
is a best fit with functional form 7(5) oc 5^'' exp(— aln^ 5). 



FIG. 3: (Color online) (Left): Log-log plot of the QFI for opti- 
mized pure input Gaussian states as a function of the noise parame- 
ter A for different values of the average photon number. From bot- 
tom to top: A'^ = {2, 5, 10}. (Right): Ratio between the Fisher 
Information of homodyne detection on coherent states and the cor- 
responding QFI, as a function of the number of photons of the 
probe states and for different values of A. From bottom to top: 
A^ = {0.5,1.0,1.5,2.0,5}. 



As we have already mentioned above, because of the form 
of the phase-diffusion map, we expect that the product ^ = 
NA plays a role in the estimation properties. In fact, by ex- 
ploring a large range of values for N and A a scaling law 
emerges from numerical analysis, which may be written as 



H{N,A) ^ k^H{N/k,kA). 



(2) 



That is, H{N,A) = N/A^{0 = N^liO/^ = ^liO/^^ 
where < 7(^) < 1 is a universal function independent 
on A and N. The larger is ^ the more accurate is the scal- 
ing law. If we fix ^ the QFI for different pairs of N and 
A have the same value, up to a rescaling by a factor fc^, 
where k is the ratio between the two average photon num- 
bers, or equivalently of the two noise parameters. The scal- 
ing is illustrated in the right panel of Fig. |2]where we report 
the quantity — ln7(^) as a function of In^ (orange points) 
together with a two-parameter fit (black curve) of the form 
liO C cxp(— a In'^ ^), that provides a good representa- 
tion of data. Using the above results, the quantum Cramer- 
Rao bound for the precision of an optimal estimator of the 
phase-shift may be written as Var(0) > ^^^^^ = ^(^^^j^i ■ 
For small values of ^ the quantity ^7(^) is of order of unity 
and thus Heisenberg limit Var(0) ~ N~'^ in precision may be 



achieved 112711 . We also found that another scaling law, though 
less accurate, holds also for the optimal squeezing fraction 



/3opt(iV,A) ~/?opt(iV/fc,fcA) 



(3) 



Though based on a physical and mathematical justification, 
we cannot expect these scaling to be exact due to the non- 
Gaussianity of the state. However they give a useful and prac- 
tical receipt to compare and predict phase estimation perfor- 
mances in different regimes of energy and noise. In the left 
panel of Fig. |3]we show the behavior of the quantum Fisher 
information at fixed average photon number as a function of 



In the noiseless case (A = 0) homodyne detection per- 
formed on input squeezed vacuum states is optimal lfl2ll . that 
is, its Fisher information F is equal to the QFI, H{N) = 
8{N'^ + N). A question thus arises on whether this results 
also holds in the presence of phase diffusion. Our numeri- 
cal findings shows that this is true for a very small amount of 
noise, i.e. A ^ 1, whereas for increasing A the ratio F/H 
is moving away from unity quite quickly. On the other hand, 
one can see that for high values of A, basically when coherent 
states are the optimal probe states maximizing the QFI, ho- 
modyne detection of the quadrature X ~ (a + at)/2 is again 
nearly optimal, i.e. its Fisher information is again approach- 
ing the value of the QFI evaluated in same conditions. In the 
right panel of Fig. [3]we plot the ratio between the Fisher in- 
formation of homodyne detection and the corresponding QFI: 
by increasing the noise A the ratio increases towards optimal- 
ity {F/H = 1). This may understood looking at the behavior 
of quadrature fluctuations AXg = (Xg) — (Xg)'^ since the 
smaller is AXg for a certain quadrature Xg, the more precise 
is the estimation of the phase-shift through this quadrature. 
In FiglH we report a contour plot of log AXg as a function 
of the squeezing fraction of the input state /? and the quadra- 
ture phase 9 for different values of A and of the overall en- 
ergy N. We see that for low noise, i.e. A <C 1, minimum 
fluctuations are obtained for the quadrature 6 — 7r/2 and 
for a squeezed vacuum state (/3 = 1), whereas after a cer- 
tain energy-dependent threshold level of noise A* = A*(A^), 
we have a jump and the minimum fluctuations are achieved 
by measuring the X quadrature (9 = 0) on coherent probes 
(/3 = 0). This behavior is different compared to the behavior 
we have obtained for the QFI, see Fig. [T] There, for inter- 
mediate values of A, the optimal squeezing fraction decreases 
monotonically from /3 = 1 to /3 = 0, whereas here we have 
only the extreme values. This exactly corresponds to the re- 
sult discussed above: homodyne detection, as far as we tune 
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accordingly the measured quadrature, is optimal for very low 
noise with squeezed vacuum probes (/3 ~ 0), and for large 
noise with coherent probes (/3 = 1), while for intermediate 
values of A homodyne detection is far from optimality. Over- 
all, we have that homodyne detection provides nearly optimal 
phase estimation in the presence of either very small or large 
phase diffusion, whereas it is still an open problem to find 
a feasible measurement attaining the ultimate precision for a 
generic value of the phase-diffusion noise parameter A. 
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FIG. 4: (Color online) Quadrature fluctuations AX| as a function 
of the squeezing fraction of the input state /3 and of the phase Q for 
different values of the noise amplitude A and the overall energy A''. 
Top left: iV = 10 and A = 0.1; top right: iV = 10 and A = 0.6; 
bottom left: A'^ = 0.1 and A = 0.1. Darker regions corresponds 
to smaller AXg. The plot in the bottom right panel illustrates the 
threshold A* (N) between the two regions where minimum fluctua- 
tions are achieved for 13 = 1,9 = tt/2 (gray area) and P = 0, 9 = 
respectively. 



In conclusion, we have attacked for the first time the prob- 
lem of finding the optimal way to estimate a phase-shift in 
the presence of phase diffusion and we have obtained the ul- 
timate quantum limits to precision for phase-shifted Gaussian 
states. By an extensive numerical analysis we have obtained 
an approximate scaling laws for both the quantum Fisher in- 
formation and the optimal squeezing fraction in terms of the 
overall total energy and the amount of noise. We also found 
that homodyne detection is a nearly optimal detection scheme 
for very small or large noise. Our results goes beyond the tra- 
ditional analysis of the quantum phase measurement problem 
and may be relevant for the development of phase-shift keyed 
optical communication schemes 1281]. 
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